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Abstract : 

A magnetic hysteresis classification of the lunar 
surface is presented. Particular emphasis is placed on 
the mode of origin of the samples which determine their 
hysteresis properties. The hysteresis ratios Rg and Rj, 
where Rjj is the ratio of remanent coercive force (%) to 
coercive force (H c ) and Rj is the ratio of saturation re- 
manence (I R ) to saturation magnetization (ls)» are the 
basic data used for the classification. 

There is a distinct correlation between natural re- 
manence (NRM) , saturation magnetization (ls)> and the hys- 
teresis ratios for the rock samples. The hysteresis clas- 
sification is able to explain some aspects of time depend- 
ent magnetization in the lunar samples and relates the 
initial susceptibility to NRM, viscous remanence, and to 
other aspects of magnetization in lunar samples. 

The framework of the hysteresis classification can 
be utilized in evaluation of the processes active in the 
formation of the regolith, and specifically in the classi- 
fication of microbreccia types. The stability of perma- 
nent remanence , the recognition of two modes of time de- 
pendent remanence, one associated with the superparamag— 
netic fraction, the other with the multidomain fraction 
in lunar samples, as well as other magnetization features, 
are compatible with the classification and indeed are 
predictable if the hysteresis properties are measured. 
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Since up to 60% of the iron in the lunar soil may be super 
paramagnetic at 4-00° K, and only 10% at 100° K, the 50% 
which becomes ferromagnetic over the cycle has the char- 
acteristics of thermoremanence and may provide for an en- 
hancement in measurable field on the dark side during a 
subsatellite magnetometer circuit. 
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Introduction; 


Permanent magnetism is present in all samples return- 
ed from the moon, and fields ranging up to 300 gamma have 
been measured by astronauts using a portable magnetometer. 
However, at present time there is no clear indication as 
to what is the magnetization mechani-sm or what is the 
source of the field which induced the permanent magnetism. 

The simplest explanation is to assume a dipole field 
in a manner analagous with thermoremanence acquisition 
in the earth's crustal layer. There are difficulties in- 
volved with assuming a dipole field. A dynamo is required 
to generate and sustain the field, and since there is no 
field at present, a dipole up to about 3 • 4x10^ years is 
necessary, active for the period of time coincident with 
the radiometric record. The apparent inhomogeneous dis- 
tribution of magnetization in the lunar surface layer 
would also seem to argue for obliteration of the dipole 
record. 

There are problems associated with remanence acqui- 
sition in iron containing samples ; The remanence acqui- 
sition efficiency of iron containing samples depends on 
the size and shape of the iron grains. Paleointensity 
estimates for a normal basalt and a contemporaneous re- 
duced sample (Wasilewski, 1972) containing the same Fe+ 
FeTi 0 ^+Fe 2 Ti 04 assemblage found in lunar rocks differed 
by two orders of magnitude. The paleointensity estimate 
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for the normal basalt was 0.33 Oersted, and for the reduced 
sample the estimate was 0.0031 Oersted. This disagreement 
was not expected even for the simple NRM-TRM comparison 
which was used for estimating the inducing field. The 
extreme case of shape was also investigated. A needle 
was given a TRM both parallel and perpendicular to the 
needle axis. For the perpendicular case no remanence 
could be measured, but for the parallel case a strong 
remanent component was measured. The only conclusion 
which can be reached is simply that the paleofield esti- 
mates published for lunar samples cannot be taken 
seriously. 

Carbonaceous chondrites and other meteoritic samples 
contain oriented vector magnetism (Wasilewski, 1972) in- 
dicating that other sources of permanent remanence induc- 
ing fields exist aside from the geomorphic dynamo. 

The body of experimental data which exists to pre- 
sent does not explain (a) the precise mode of remanence, 

(b) the type of inducing field, (c) the time of magneti- 
zation, (d) the magnitude of the inducing field. To at- 
tempt to answer the numerous questions precipitated by 
the measurement of permanent remanence in the lunar sam- 
ples, we must first understand what kind of ferromagnetic 
dispersion we are dealing with and what effect thermophy- 
sical processes active in a lunar environment would have 
on the remanence characteristics of lunar samples. 
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Magnetic hysteresis characterization provides a step 
in the proper direction. In this paper I will demonstrate 
that the magnetic hysteresis characteristics are diagnos- 
tic and can he utilized as a classification tool for all 
lunar samples and that thermo physical processes will ef- 
fect changes in the hysteresis characteristics of lunar 
samples which are both predictable and useful. 

Some Magnetic Hysteresis Aspects of Iron Dispersions Rele- 
vant to Lunar Research : 

The makeup of the ferromagnetic dispersion in terms 
of the size and shape of the ferromagnetic components and 
the range of saturation magnetization will determine the 
magnetic hysteresis of the dispersion. 

If a superparamagnetic dispersion of iron has a mag- 
netization value, M s , at saturation, it will after time, 
t, have a value, M=M r , defined by 

M r =M s EXP(-t/T ) 1 

where l/f = f 0 EXP(-KV/kT) (Neel, 1949). A spherical 
iron particle of radius, 115 &, with the crystalline ani- 
sotropy, Ki, being the only anisotropy, has a value of 

]_0“1 seconds, and when the radius increases to 150 X, 
the value of T = 10 9 seconds, and the particle is a stable 
ferromagnetic grain. Superparamagnetism is defined with 
H c =0, and though it does not contribute to H c , it will 
reduce H c , if present with single domain particles, ac- 
cording to the argument H c *q/lR in the equation 


- 5 - 



-1 


2 


H C =H C l+(H c .q/l R ) €/(l-€) 
where q=V M x /3kT for superparamagnetism, where V is the 
volume, M fl , the magnetization, k, Boltzman's constant, 
and T, the absolute temperature. If large multidomain 
grains are present with single domain grains, the same 
argument holds as above, except that q=(NM s ) -1 , where N 
is the demagnetization factor. The above discussion is 
based on the work of Kneller and luborsky (1964), sum- 
marized by Kneller (1969). 

It has been experimentally verified by Wasilewski 
(1972 ) that mixing two distributions represented by 
discs of basaltic rock, one with H c »50 Oersted, the 
other with H c =420 Oersted, produces a reduced H c value 
in accordance with the calculations of Wohlfarth (1954) 
for two sorts of particles with H Cl <H C2 . Exact agree- 
ment cannot be expected since each H c disc value is the 
effective H c value for a distribution itself, with the 
resultant distribution being bimodal. 

In a fine dispersion, the initial susceptibility 
is more sensitive to the larger particles and the approach 
to saturation is largely governed by the smaller particles 
in the distribution. 

Since a single domain particle cannot change its 
magnetization except by rotation of the spontaneous 
magnetization against anistropy forces, the particle 
has coercivity determined by the crystalline anisotropy 

characteristic of the material. For iron and dilute 
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FeNi alloys with high saturation magnetization the co- 
ercivity is proportional to the magnetization, and it 
is easily demonstrated experimentallly and theoretically 
(Kittel, 1949 ) that shape anisotropy is the usual ani- 
sotropy that figures in the evaluation of the coercivity 
of irregular grain dispersions. 

For oriented iron particles (Luborsky and Paine, 1962) 
the H Q value and the Rj value depend on the measurement 
angle, 9, with respect to orientation. The maximum Rj 
occurs for 9=0° and the minimum for 9=9Q°. The maximum 
H c value occurs at an angle, 9=50°, for dilute assemblages, 
but for increased packing magnetostatic interaction becomes 
important and the 50° maximum is suppressed. The minimum 
values of H c occur at 90°. The lunar soil may contain 
an oriented thin film surface of metallic iron, and it 
is clearly shown (Carter, 1971) that the surface of the 
soil particles contain fine (0.03 to 0.5 /*m) mounds of 
magnetostatically interacting iron. A definite surface 
anisotropy which should provide a magnetic stratigraphic 
boundary zone in a core should be present, particularly 
if the metallic coating is associated with discrete 
cratering events. 

Rotational hysteresis (W r ) measurements provide 
a useful check on the existence of exchange anisotropy 
and magnetically uniaxial components. The only rotation- 
al hysteresis work to date was performed by Runcorn et al. 
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(1970). They wrongly ascribe the large W r values at high 
measuring fields to iron, explaining the W r losses to 
shape anisotropy. Their Figure 9 appears to be impossible 
as the work of Jacobs and Luborsky (1957) and Campbell 
et al. (1957) clearly demonstrates. High field W r is due 
to anisotropy of uniaxial, induced, or exchange character 
(see Meikel john, 1962). For Ni, Fe^O^., and Fe in bulk 
or dispersed systems, no significant high field W r is 
expected. That Runcorn et al. have observed high field 
W is of course significant, but their interpretation is 
quite wrong. They also observed high field W r for an as- 
semblage of polycrystalline iron spheres (5p m) . There 
is something peculiar with their experiment since their 
result (a) does not agree with theoretical studies nor 
(b) does it agree with previously determined W r -H curves 
for single and multidomain iron. Intergrowths such as 
Fe-FeS, where exchange anisotropy has been observed 
(Greiner et al. , 1961), or Fe-FeO epitaxial overgrowths 
will give the results they observed. This proposition 
can be easily verified (Meikeljohn, 1962). 

In iron dispersions, if size is important, then both 
H c and Hr are linear functions of temperature. If shape 
is important, there should be different effects, since 
the temperature coefficient of H c depends on shape. 

There are two types of time dependent effects which 
can operate in the lunar samples. The first is ascribed 
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to the superparamagnetic component which was considered 
earlier (equation 1 ), and the other is due to multi- 

domain iron which is described by Al=S log t+constant. 

It should be noted that S, the viscosity coefficient 
depends on (a) the shape of iron particles (Street e_t al. . 
1952; and Yakubaylik, 1967) and on (b) the alloy composi- 
tion, since the magnetization value and structure will 
change . 

Magnetic Hysteresis in Lunar Samples ; 

A tabulation of some data presented by Nagata £t al . 
(1970, 1971, 1972) for Apollo 11, 12, and 14 samples is 
presented (Table 1). The samples are placed into two 
groups based on Rj and Rjj values. Group A contains crystal 
line rocks and thermally metamorphosed breccia samples, 

(P 4 - Jackson and Wilshire, 1972). Group B contains 
the lunar fines and the welded breccia samples, (P 1 - 
Jackson and Wilshire, 1972). 

The hysteresis loop and the various parameter^ includ- 
ing the measurement paths, are illustrated in Figure 1. 

The hysteresis ratios Rj and R H are demonstrated to be 
characteristic for each group of natural materials 
(Wasilewski, 1970, 1972 , 1972 , 1972 ). The Rj and 
R h values plot as two distinct groups (Figure 2k). One 
group with Rj<0.02 and R H <10.0 contains the igneous 
rocks and thermally metamorphosed breccia samples, the 
other group with Rj>0.04 and Rh> 10.0 contains the fines 
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and welded breccia samples. 

The H c values for all samples fall between 10 and 
50 Oersted (Figure 2B), as shown by the dotted lines, and 
the lunar samples are subdivided into two groups on the 
basis of the Rj vs. H Q data, one group with Rj between 
0.004 and 0.02 contains the crystalline rocks and thermal- 
ly metamorphosed breccia samples, the other with Rj between 
0.04 and 0.072 contains the lunar fines and welded breccias. 
The Rj vs. H c data clearly demonstrate that H c isa meaning- 
less parameter when applied to lunar samples, and the rea- 
sons for this are quite obvious. Wasilewski (1972 , 1972 ) 
has discussed this aspect of natural materials in terms 
of the experimental and theoretical work of Meikeljohn 
(1953), Wohlfarth (1958), Stoner and Wohlfarth (1948), 
and Kneller and Luborsky (1964). Mixtures of single domain 
iron particles with either multidomain or super paramagne- 
tic particles will result in a reduction in H c . A given 
volume fraction of super paramagnetic material will reduce 
H c more than an equivalent volume of multidomain material. 

The Rj vs. Hr plot (Figure 2B) results in two dis- 
tinct and characteristic groups defined on the basis of 
Rj and Hr values. For the first group, which includes 
the crystalline rocks and thermally metamorphosed breccia 
samples, Rj<0.02, and Hr varies between 75 and 180 Oersted, 
while for the second group, which includes the fines and 
welded breccias, Rj>0.04, and Hr varies between 300 and 
520 Oersted. The Rj vs. Hr plot appears to be most infor- 
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mative 


The presence of three discrete iron modes provides 
for a lCr range in grain H c ranging from H c =0 for the 
superparamagnetic material to Hp^lO for single domain 
material. All of the the hysteresis loops thus far pre- 
sented for lunar samples show constriction in the region 
of low measuring fields. This constriction can he found 
in partially oxidized terrestrial basalts (Wasilewski, 
1972), chondrite meteorites (Wasilewski, 1972), and in 
the lunar samples (Nagata et al . , 1970, 1971, 1972). 
Discussion : 

The data presented in Table 1 is quite systemmatic 
despite the limited sampling. The samples were arranged 
in order of increasing Rjj values, and all other parameters 
follow a systemmatic increase or decrease accordingly. 

It is also of interest to note that Nagata has defined 
two types of viscous (VRM) behavior for the analyzed lunar 
samples, and type I VRM behavior is associated with group 
A hysteresis characteristics, and type II VRM behavior 
is associated with group B hysteresis behavior. 

Saturation in the positive H direction will result 
in a finite remanence when H is returned to zero. This 
remanence depends on the magnetization mechanisms in the 
single domain and multidomain size fractions, but does 
not depend on the superparamagnetic or paramagnetic com- 
ponents. The H c value is meaningless, as was indicated 
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earlier, and the H R field is that required to produce a 
zero remanence state in the sample after saturation. A 
direct sample to sample comparison of this value is in- 
formative since it indicates the magnitude of the field 
needed to reverse half of the irreversible remanence 
left in the sample after saturation. The Rj value is al- 
so a relative indication of the amount of irreversible 
magnetization, and, as can be seen in the Rj vs. 
plot, this value increases as H R increases. 

In practice the R^ value is a complex number, in that 
it can depend on 

(a) grain size 

(b) grain shape 

(c) packing fraction 

(d) grain alignment 

(e) exchange anisotropy 

(f) degree of strain anisotropy 

(g) degree of strain anisotropy 

(h) anisotropic structures 

(i) size distribution 

(3) amount of superparamagnetism. 

The lower the H R value the more viscous is the size 
fraction larger than single domain size. There are two 
components which are included in time dependent changes, 
(a) the superparamagnetic fraction and (b) the multidomain 
fraction. This is reflected in the Rjj values, the H c 



values, and the initial susceptibility. For a constant 
Hr value the Rr value will depend on the amount of super- 
paramagnetism. A more detailed study of viscous behavior 
based on low temperature analysis of the blocking spec- 
trum will allow the size spectrum of superparamagnetism 
to be separated from the multidomain size spectrum. 

This is a critical aspect of lunar magnetization 
since the total magnetization of the sample depends 
on the handling of the sample over a temperature range 
equivalent to the range in which the Apollo samples 
were handled, and unless the blocking spectrum for 
this range is defined I do not believe we can eliminate 
the noise associated with the transfer of the samples 

I 

from the moon to earth. 

Based on the experimental studies of Nagata et al . 
(1970, 1971, 1972), Thorpe and Sentfle £t al. (1970, 
1971), Jedwab (1971, 1972), Carter (1971), Carter and 
MacGregor (1970), Chao at al. (1971), and many others 
we can conclude that the relative significance of super- 
magnetism (SP) can be specified: 

SP(Fines)> SP(Welded Breccia)> SP(lgneous Rocks) >' 
^ SP(Thermally Metamorphosed Breccia). 
Therefore, the H c value will be reduced relatively more 
significantly in the reverse order of the list above, 
and the Rr value will depend on the Hr value almost 
exclusively, particularly since the H c values (Table 1 
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and Figure 2B) are limited to a narrow range. 

The H r value will depend on the shape and range of 
the size distribution. For igneous rocks and thermally 
metamorphosed breccia samples the H R value is expected 
to be lowest since the majority of iron grains will be 
multidomain with single domain or small multidomain 
grains carrying any stable component. As the shape of 
the size distribution becomes wider, shifting to a smal- 
ler size range, or bimodal, as in the case of the fines 
and welded breccia, the H R value should increase since 
a large fraction of the iron will be single domain or 
nearly so. The fines contain the broadest size distri- 
bution of ferromagnetic components, and welding the fines 
should not produce any significant alterations to the 
distribution. Thermal metamorphism will destroy some 
or all of the superparamagnetism, and an effective shift 
of the size distribution will result in a decrease in Hr, 
R R ,and Rj, as observed. 

The data thus far accumulated provides distinct 
evidence that the lunar samples are quite different from 
terrestrial samples, mainly because there is iron in the 
lunar samples. 

The consistent R R and H R values for the lunar fines 
suggest a common origin and common size and shape modes 
for the ferromagnetic fraction. Though data is sparse 
at present, there is a definite suggestion that the per- 
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manent remanence of each lunar sample can he understood 
in terms of the hysteresis properties of each sample. 

The path the hysteresis parameters will take depends on 
the mode of origin and acting thermomechanical processes. 

A soil can he welded, thermally metamorphosed, or shock 
lithified. Each of the three processes will produce 
different effects depending on 

(a) the behavior of the SP component 
(h) the change in the ferromagnetic component size 
distribution 

(c) the magnetostatic interaction 

(d) the influence of the shape of the ferromagnetic 
components 

(e) the composition of the alloy, i.e. Ni+Co content. 
Conclusions : 

I. There are well defined groups of lunar materials 
based on their magnetic hysteresis properties. 

II. Natural remanence correlates directly with satura- 
tion magnetization, the Rj value, and the initial 
susceptibility for crystalline rocks. 

III. In lunar samples the low H c values are due in part 
to the presence of superparamagnetic iron and in 
part to multidomain iron. 

IV. The H c value is a meaningless parameter for lunar 
samples. 

V. The H r value is important and diagnostic since it 
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is due to near single domain and multidomain grains. 
It is independent of the super paramagnetic fraction. 
YI. Microbreccia samples can be analyzed non-destructive- 
ly and classified as to their mode of origin based 
on magnetic hysteresis. This approach will be more 
effective with a more systemmatic array of data. 

VII. It is possible that the measured surface fields 

can be due to a magnetized regolith. This requires 
further verificaU on, but it is a distinct possibil- 
ity which must be considered. 

VIII. A soil can be subjected to welding, thermal meta- 
morphism, or shock lithification. Each of the three 
processes will produce different hysteresis charac- 
teristics and remanent states depending on 

(a) the behavior of the SP component 

(b) the change in the effective grain size distri- 
bution 

(c) the development of eutectic intergrowths 

(d) the phase composition, i.e. Ni+Co content 

(e) the thermomagnetic anomalies of the ferro- 
magnetic components. 

IX. Two distinct viscous or time dependent features 
are characteristic of iron in the lunar samples. 

The first concerns the superparamagnetic fraction, 
and the second concerns the soft multidomain iron 
fraction. Both features are temperature dependent. 
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X. Prom the hysteresis data and other published 
results it appears that relative superpara- 
magnetism (SP) can be specified: 

SP(Pines)> SP(Welded Breccia)^ SP( Igneous Rocks) 
^ SP(Thermally Metamorphosed Breccia). 

The H c value will be relatively reduced in the 
opposite order. 
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TABLE 1 

(Tabulation of data from Nagata et al. 
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FIGURE CAPTIONS 


Figure 1 - Schematic diagram of magnetic hysteresis loop 

illustrating the measurement paths and defining 
the loop parameters 

Ir - Remanent magnetization 

Ig - Saturation magnetization 

H c - Coercive force 

Hr - Remanent coercive force 

R i - Ir/Is 

% - %/H 0 

(1-2-1, measurement path for full loop; 1-3-0, 
measurement path which reduces Ir to zero and 
defines Hr) 


Figure 2A - Relationship between Rt and Rr for lunar samples 
(Data from Nagata et at., 1970, 1971, 1972) 

Figure 2B - Relationship between Rt and Hr for lunar samples 
(Data from Nagata £t at., 1970, 1971, 1972) 

The Rt vs. Hq values Tor the lunar samples are 
included in the dotted areas labelled I and II. 
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